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We perform a detailed analysis of above threshold ionization (ATI) in atoms within the strong 
field approximation (SFA) by considering spatially inhomogeneous monochromatic laser fields. The 
locally enhanced field induced by resonance plasmons is an example for such inhomogeneous fields. 
We investigate how the individual pairs of quantum orbits contribute to the photoelectron spectra 
and the angular electron momentum distributions. We demonstrate that the quantum orbits have 
a very different behavior in the spatially inhomogeneous field when compared to the homogeneous 
field. In the case of inhomogeneous fields, the ionization and rescattering times differ between 
neighboring cycles, despite the field being monochromatic. Indeed, the contributions from one cycle 
may lead to a lower cutoff, while another may develop a higher cutoff. Within our model, we show 
that the ATI cutoff extends far beyond the semiclassical cutoff, as a function of inhomogeneity 
strength. Furthermore, the angular momentum distributions have very different features compared 
to the homogeneous case. For the neighboring cycles, the electron momentum distributions do not 
share the same absolute momentum, and they do not have the same yield. 
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I. INTRODUCTION 

In the context of the interaction of matter with strong 
laser fields, above threshold ionization (ATI) [l], [2[ has 
attracted considerable interest, since the behavior of the 
laser-ionized electrons during ATI serve as a very valu- 
able tool for laser pulse description. The electric field in a 
few-cycle pulse can be characterized by its duration and 
by the so-called Carrier Envelope Phase (CEP) [j| |j]. 
In particular, for a better control of the system on an 
attosecond time-scale it is important to find reliable and 
direct schemes to measure the absolute phase of few-cycle 
pulses. The investigation of ATI generated by few-cycle 
driving laser pulses plays a key role in the CEP charac- 
terization due to the sensitivity of the energy and angle- 
resolved photoelectron spectra to the value of the laser 
electric field absolute phase [5|, Q . In addition, it seems 
that the high energy region of the photoelectron spectra 
is most sensitive to the absolute CEP and consequently 
ATI yields with larger cutoffs is needed in order to de- 
scribe it [HE}- 

Recent studies have demonstrated that the cutoff of 
the high harmonic generation and photoelectron spec- 
tra could be extended further by employing spatially 
inhomogeneous fields [a-tla]. In particular, recent ex- 
periments using a combination of rare gases and plas- 
monic nanostructures have demonstrated that the har- 
monic cutoff of the gases could be extended further by 
using the locally enhanced field induced by coupling laser 
pulses to nanosvstems [lil Il5| . In such system, the locally 
enhanced field is not spatially homogeneous in the region 
of electron dynamics. This is a result of the strong con- 
finement of the plasmonics spots and the distortion of 
the electric field by the surface plasmons induced by the 
nanosystem. In addition, recently, solid state nanostruc- 



tures have been employed as a target to study t he p ho- 
toelectron emission by few intense laser pulse [13, [l8| . In 
this process, the emitted electrons have ener gy far be- 
yond the usual cutoff for noble gases (see e.g. |18l423j ). 

Hitherto, the theoretical approaches, including solv- 
ing the Time Dependent Schrddinger Equation (TDSE) 
and the Strong Field Approximation (SFA), for studying 
strong field phenomena are mainly based on the assump- 
tion that the laser electric field is spatially homogeneous 
in the region of electron dynamics [2J, (25[. This, hy- 
pothesis, however, is no longer valid for the locally en- 
hanced plasmonic field. Indeed, in such system the strong 
confinement of the electrons in the plasmonic hot spots 
generates a spatially inhomogeneous inhomogeneous elec- 
tric field. It means the driven electron will experience a 
Lorentz force which depends on position and its subse- 
quent motion will be strongly affected. As a result, new 
physics will emerge in the interaction between matter and 
strong laser fields. Until now, few theoretical approaches 
have been developed investigating the strong field phe- 
nomena in such spatially inhomogeneous fields |10Hl3j. 

In the present paper, we employed SFA to investigate 
the ATI by resonant plasmon field enhancement. In ad- 
dition, we use a SFA based on saddle point methods 
rather than a full numerical SFA approach to obtain the 
ATI photoelectron spectra and electron momentum dis- 
tributions. The saddle point methods lead to equations 
that can be directly related to the classical equations 
of motion of an electron in a laser field. As a result, 
they provide a space-time picture which gives us addi- 
tional physical insight. In this work, which based on 
this method, we investigate the individual electron tra- 
jectories and demonstrate their contributions to the ATI 
spectra and the electron momentum distributions. In 
addition, since the imaginary part of the saddle point 



equations can be related to the width of the potential 
barrier through which the electron tunnels, we examine 
the ionization probability of the electron for each trajec- 
tory. 

This article is organized as follows. In Sec. [TT1 we 
present the SFA transition amplitudes for both direct 
and rescattered ATI driving by lineally spatially inhomo- 
geneous field. We start from common expressions based 
on homogeneous fieldsf Sec. Ill A|) . and subsequently show 
the modified transition amplitude by considering inho- 
mogeneous fields with linear dependence (Sec. Ill B[) . In 
Sec. 1111 A[ we discuss the saddle-point equations and an- 
alyze them in terms of quantum orbits in parallel to 
their classical counterpart electron trajectories, in terms 
of both electron energy and momentum . In the next sec- 
tion, Sec lIIIBl we present the ATI spectra based on the 
analysis given in Sec MI A T] In Sec. lIIICl we provide an- 
gular momentum distributions of the direct and rescat- 
tered ATI based on the analysis given in Sec. IIII A 21 
Finally, in Sec. IIVI we conclude the paper with a few 
summarizing remarks. 



II. THEORY 



A. Transition Amplitude 



Thereby, k, p and I p denote the intermediate and final 
momentum of the electron and the ionization potential 
of the field-free bound state \<po), respectively 

The prefactors of the transition amplitudes Mdir and 
M res read 
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and 



Vp k = (p(t)\V\k(t)) 

= ^p/ d3rex PH(p-k)-r]y(r)0 o (r) (7) 

where V gives the interaction of the system with the laser 
field. Here, k(t') = k and p(t') = p in the velocity 
gauge, and k(t') = k + A(t,) and p(t') = p + A(t') in 
the length gauge. 



The strong field approximation (SFA) is based on two 
assumptions, namely (i) the influence of the laser field 
is neglected when the electrons are bound to their tar- 
get atoms and (ii) the binding ionic potential is neglected 
when the electrons are in the continuum. As a result, the 
free electrons in the continuum are described by field- 
dressed plane waves, known as Volkov states [26|, I27J . 
Based on the assumption that the laser electric field 
does not change with respect of the position in the re- 
gion where the electron dynamic takes place, the well- 
established Keldysh- Faisal- Reiss (KFR) model [H-[3lj], 
gives the SFA transition amplitude for the direct and 
rescattered ATI as (in atomic units): 
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and 



M rea = - dt dt 1 d 3 kV pk V w e lSo ^^'\ (2) 
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where the corresponding actions are given by 

So(f) = - J X [V + f T)? dr + I p t (3) 



and 



(4) 



B. Strong Field and Newton's Equations 

We now study a case in which the laser field has a spa- 
tially inhomogeneous character in the region where the 
electron motion takes place. We start by examining how 
the action of the SFA is connected to the classical elec- 
tron trajectories. In the SFA, the action is defined in 
terms of the vector potential field, which is the counter- 
part of the velocity x(t) in Newton's equation of motion. 
The laser potential Vl due to the laser field E(t, x) is 
defined as 



V L =xE(t,x), 



(8) 



and the Newtonian equation of motion for an electron in 
this field is given by 



x(t) = V X V L = -xV x E{t, x) - E(t, x). 



(9) 



For the homogeneous case, in which the laser field does 
not have spatially dependency and the laser electric field 
is just E(t), the Newton equation reads x(t) = —E{t). 
On the other hand, if the spatial dependence of the en- 
hanced laser electric field is perturbative and linear with 
respect to position, then the field can be approximated 
as 



E(t,x) ~E(t)(l + ex), 



(10) 



where e <C 1 is a parameter that characterize the strength 
of the inhomogeneity. 



Indeed, the above approximation corresponds to the 
first term of the actual field of a plasmonic nanostructure 
with spherical shape (321 • By substituting Eq. (flO|l into 
©, we have 



x(t) = -E(t)(l + 2ex{t)). 



(11) 



This is the effective laser electric field that the electron 
feels along the trajectory x(t), which describes its motion 
once it is laser freed to the continuum. 

Classically, the electron trajectory can be found by 
solving Eq. (|lip . In here, we solve it by applying the 
Picard iteration 33[ method and restrict ourselves to the 
first order (for more details see [12[ ). Based on the con- 
dition that the electron starts its movement at the origin 
with zero velocity, i.e. x(0) = and x(0) = v(Q) = 0, we 
obtain 



x{t) = f3{t) - p{t ) - A{t ){t - t ). 



(12) 



with f3(t) = J Q dt'A(t'). In addition we assume that at 
time to the potential field is zero, then 



x(t) = / dt'A(t'). 



(13) 



By using Eqs. (TTTjl and (JT5]) the effective vector poten- 
tial along the electron trajectory A tr (t, x) reads 



where 



and 



A tr (t)=A(t)+2eA c (t), 



t j-t 



(14) 



A c (t) = / dt"A{t") - / dt"A 2 {t"), (15) 



A{t) = - / dt'E(t'). 



(16) 



The next step is to incorporate the above defined in- 
homogeneous vector potential field into the general ex- 
pression of the transition amplitudes, i.e. Eqs. (fTJ) and 
(|2|) for the direct and rescattered ATI, respectively. Con- 
sequently, we replace the vector potential A(i) in M^ir 
and M res by the one defined in Eq. (fT4"|) . respectively. As 
a result, the modified actions yield 



S(t') = S Q (t')-2e I A c (r)[p + A(r)]dr 



, Al{r)dr 



(17) 



and 



A c (r)[p + A(r)]dr 

-2e M c (r)[k + A(r)]dr-2e 2 / A 2 c ( T )dT (18) 

where So(t') and 5o(k, t,t') are defined in Eqs. © and 
Q, respectively 



C. Saddle-point equations 

For large driving-field intensities, the transition ampli- 
tudes Eqs. (fTJ) and @ are strongly oscillatory integrals 
which can be evaluated using the saddle-point or steep- 
est descent method [3J, |35[ . The solutions of the saddle 
point equations are directly related to the classical tra- 
jectories, which allow us to compare them with the quan- 
tum orbits. This method requires obtaining the saddle 
points where the action Eq. ([3]) for the direct ATI and 
the action Eq. (jj) for the rescattered ATI are station- 
ary, i.e. for which d t S(t') = d t S(k,t,t') = d t ,S{k,t,t') = 
dkS(k,t,t') = 0. In this paper, we use a specified steep- 
est descent method called uniform approximation to take 
care of those saddles points which are not well separated 
(for a detailed discussion see Ref. [3^]). For the direct 
ATI, the stationary condition upon t' lead to the saddle- 
point equation 



[p + A(i')] +4eA(0 



-2L, 



with 



X(t')=A c (t')[p + A(t')] + eA 2 c (t') 



(19) 



(20) 



For the rescattered ATI, the stationary conditions 
upon t, t' and k lead to the saddle-point equations 

[k + A(t')] 2 + 4e«(t') = -2I P (21) 



/ dr[k + A(t)] + 4er/(r) = (22) 

Jt' 



and 



[p + A(i)] 2 _ [k + A(t)] 



2eA c (i)[k-p] (23) 



with 



and 



K(t')=A c (t')[k + A(t')]+eA 2 c (t'), 



tj(t) = / A c ( T )d7 



(24) 



(25) 



Equations (|2Tj) and (fT9|) express the conservation law of 
energy for the electron tunnel ionized at the time t' . Fur- 
thermore, Eq. (|22p constrains the intermediate momen- 
tum of the electron and guarantees that the electron re- 
turns to its parent ion. Finally, Eq. ((23)) gives the energy 
conservation of the electron at the time t, when upon its 
return inelastically rescatters with the core. 

The inhomogeneous character of the laser field gives 
the additional terms A(i'), n(t') and r](t) in the Eqs. (|19[) . 
([2~T1) . and (|22l) . respectively. These terms vanish for the 
homogeneous case, i.e when e = 0. Furthermore, for the 
homogeneous case Eq. (f2"Tj) has no real solutions, unless 



Ip — > 0, due to tunnel ionization, which has no classi- 
cal counterpart. Therefore, the solutions of the saddle 
equations are generally complex. For the inhomogeneous 
case, however, it is not very upfront to constrain the lim- 
its, in which the solutions of Eq. (f2~Tj) are real. Neverthe- 
less, in our case, the electron will most likely reach the 
continuum via tunnel ionization, since e is a very small 
parameter. It means that the solutions of these saddle 
point equations are still expected to be complex. In addi- 
tion, the maximum kinetic energy that the electron gains 
in the continuum now depends on the nonhomogeneous 
character of the field (i.e. on the value of e) and A c (t), 
and it does not have anymore the conventional value of 
3.17£/n, where U p = Eq/(4u; 2 ) is the ponderomotive en- 
ergy m M- 



III. RESULTS 
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A. Quantum orbits 

In this section, we study the solutions of the saddle- 
point equations to examine the role of individual trajecto- 
ries in the ATI photoelectron spectra cutoff and angular 
momentum distributions by performing a quantum-orbit 
analysis of the problem. The concept of the quantum- 
orbits emerges from the fact that the solutions of the 
saddle-point equations are related to the classical trajec- 
tories of the electron in the laser field and, in addition, 
provides information on quantum aspects such as tun- 
neling and interference. To get a better insight into the 
inhomogeneous case, we employ a monochromatic field: 



E(t) = E sin(wt)f 



(26) 



where e z is the polarization vector along the z-axis (see 
Fig.CE]). 

By using the relationship defined in Eq. (fT4|) and ap- 
plying some trigonometric identities, the laser effective 
vector potential field associated with the monochromatic 
field Eq. (|2"6")l along the electron trajectory reads 



A tr (t) = A cos(wt) + 2eAc(t), 
where Aq — Eq/lo and 

A c (t) = A 2 Q sin(cji)/4w - A%t/2. 



(27) 



(28) 

If we consider the limit I p — > 0, i.e. the electron 
reaches the continuum with zero kinematical momentum, 
the drift momentum of Eq. (|2 1|) written in terms of the 
ponderomotive energy yields 



k = -2^/C/p cos(wt') + e — sin(wi') - 2U p t' . (29) 



1. In terms of electron energy 

We consider the final momentum p is parallel to the 
laser field and write it in terms of kinetic energy to illus- 
trate the simple-man model for ATI. We solve the saddle 



FIG. 1: (color online) Schematic representation of the laser 
electric field E(t) and the corresponding vector potential field 
A(£) for a monochromatic field defined by E(£) = Eo sin(u)t)e x 
(we show only 3 cycles) vs. t in a.u. The arrow indicates 
the classical times around which the electrons leave to the 
continuum and return to the core (approximately the field 
maxima and the field crossings, respectively). The pairs of the 
orbits are indicated by the labels Pair n , where n ranges from 
1 to 4. The fields are normalized to A(t)/Ao and E(t)/i?o , 
where Aq and -Eo are the field peak amplitudes. 



point equations defined in Eqs. (TT9")) - (|23[) in terms of the 
ionization t' and rescattering t times, respectively. For 
more close analysis, we restrict ourselves just to the so- 
lutions of the first 3 cycles of our defined monochromatic 
field, as shown in Fig. [TJ Classically, it is most probable 
that the electron ionizes at the electric field maxima and 
returns to its parents ion at the electric field crossings, 
i.e. when E(t) — 0. In Fig. [lja) we depict the cases 
when the electron returns to the core at a time about 7r 
(t = it/uj sa 55 in a.u) while it tunnels at the field max- 
ima at times about — 7r/2 (t sw —28 a.u.), —3tt/2 (t ra —83 
a.u.), -5tt/2 (t « -138 a.u.), and -7tt/2 (t w -193 a.u.). 
A similar analysis was carried out in Ref. 39] for an ho- 
mogeneous field. Figure QJb) shows the cases when the 
electron returns to the core at a time about — -n (t « —55 
a.u.) while it tunnels at the field maxima at times about 
tt/2 (t « 28 a.u.), 37r/2 (t w 83 a.u.), 5tt/2 (t w 138 a.u.), 
and 7tt/2 (t w 193 a.u.). 

In Fig. [21 we plot the energy of the rescattered electron 
as a function of the real part of the ionization t' and 
rescattering t times for the case with e = (panels (a) 
and (b)), e = 0.003 (panels (c) and (d)) and e = 0.005 
(panels (e) and (f)). The panels in the left and right 
columns show the solutions correspond to the pairs given 
in Fig.QJa) and Fig.QJb), respectively. 

For the homogeneous case, the ionization and rescat- 
tering times corresponding to each cycle are identical, i.e. 
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FIG. 2: (color online) Energy of the rescattered electron as a 
function of the real part of the release time t' and the rescat- 
tered time t of the electron for all given pairs in Fig. JT}. We 
consider an hydrogen atom, for which the ground-state en- 
ergy is I p = —0.5 a.u., in a linearly polarized, monochromatic 
field of frequency uj = 0.057 a.u. (corresponding to a wave- 
length A = 800) and intensity I = 3 X 10 14 W/cm 2 . Panels 
(a) and (b) correspond to e = 0, while panels (c) and (d) 
show e — 0.003 and panels (e) and (f) depict e — 0.005. The 
left and right columns show the solutions corresponding to 
the pairs given in Fig. QJa) and Fig. \Hb), respectively. The 
dashed and solid lines correspond to the long and the short 
orbits, respectively. 



for the pairs when electron leaves at about — (n + 1/2)tt 
and returns at about n (panel (a)) and for the pairs when 
it leaves at about (n + 1/2)tt and returns at about — tt 
(panel (b)). However, this is not true for the spatial 
inhomogeneous case. In here, for the case when elec- 
tron leaves at about — (n + l/2)7r and returns at about 
tt (panels (c) and (e)) the energy cutoff is larger, while 
for the case when electron leaves at about (n + 1/2)tt 
and returns at about — tt (panels (d) and (f)) the energy 
cutoff is much smaller. In general, as a function of inho- 
mogeneity of the field e, for the former case the energy 
cutoff will be extended, while for the latter case the en- 
ergy cutoff will move toward lower energy values. For 
the homogeneous case, the shortest pair (Pairi) has the 
largest cutoff, while for inhomogeneous case, the e will 
determine the pair with largest cutoff. For instance, for 
e = 0.005 (panel (a)) Pair^ has larger cutoff than Pain. 

Figure |3] depicts the energy of the rescattered electron 
as a function of the imaginary part of the release time 
t' of the electron, for the case with e = (panels (a) 
and (d)), e = 0.003 (panels (b) and (e)) and e = 0.005 
(panels (c) and (f)). The upper and lower panels show 
the solutions correspond to the pairs given in Fig. Q^a) 
and Fig. [TJb), respectively. 

For the homogeneous case, the imaginary part of the 
release time t' corresponding to each cycle are identical, 



12 

9 

6 

3 


12 

9 

6 



0l — 
0.5 




FIG. 3: (color online) Energy of the rescattered electron as 
a function of the imaginary part of the release time t' of the 
electron for the same parameters as in Fig. [2] and for all given 
pairs of Fig. [I] Panels (a) and (d) correspond to e = 0, while 
panels (b) and (e) show e = 0.003 and panels (c) and (f) de- 
pict e = 0.005. The upper and lower panels show the solutions 
correspond to the pairs given in Fig. [Ha) and Fig. [T^b) , re- 
spectively. The dashed and solid lines correspond to the long 
and the short orbits, respectively. 



i.e. for the pairs when electron leaves at about —(n + 
1/2)7T and returns at about 7r (panel(a)) and for the pairs 
when it leaves at about (n + l/2)7r and returns at about 
—tt (panel(d)). However, this is not anymore valid for 
the inhomogeneous cases. The inhomogeneity strength 
of the field e strongly affects the imaginary part of the 
release time t', and as a result, the ionization rate of the 
electron. 
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FIG. 4: (color online) Energy of the direct electron as a func- 
tion of the real and imaginary part of the release time t' for 
the same parameters as in Fig. [2] Panels (a) and (d) corre- 
spond to e = 0, while panels (b) and (e) show e = 0.003 and 
panels (c) and (f) depict e = 0.005. The upper and lower 
panels show the real and imaginary time of the ionization, re- 
spectively. The black and red lines show when electron leaves 
the atom at 7r/2 and 37r/2, respectively. 

In Fig. [U we present the energy of the direct electron 
as a function of the real and imaginary part of the release 



time t'. For the homogeneous case, the imaginary part 
of the ionization time is the same for the case when elec- 
tron leaves at ir/2 and 3tt/2 as shown in panel (d). For 
this case the real time of trajectories are symmetric at 7r 
(panel (d)). On the other hand, for the inhomogeneous 
case, the imaginary time of the ionization is different for 
the case when electron leaves at ir/2 and 3tt/2 as shown 
in panels (e) and (f). The deviation becomes larger as 
function of e. In addition, the real part of the saddle 
point equations when the electron tunnels at tt/2 and 
3ir/2 become less symmetric at around tt (panels (b) and 
(c)). 



2. In terms of electron momentum, 




2 2 4 6 

P|i /[u p r 



We now investigate the solutions of the saddle point 
equations defined in Eqs. (|19[1 - (|23[) in momentum space. 
In order to analyze this, we rewrite Eqs. (fT9f and (|23|) in 
terms of the electron momentum components parallel pn 
and perpendicular p± to the laser-field polarization, 



[PH+A(0] +4e\(t') = -2I p -p 2 ± -4eA c (t')p\\ (30) and int(msity , =: 
with 



FIG. 5: (color online) Real part of the tunneling t' and rescat- 
tering t times for the electron, as a function of its parallel 
momentum p« , when the electron tunnels around tt/2 and 
rescatters at about 2tt. We consider an hydrogen atom, for 
which the ground-state energy is I v = 0.5 a.u., in a linearly 
polarized, monochromatic field of frequency ui = 0.057 a.u. 



3 x 10 14 W/cm 2 . Panels (a) and (d) cor- 



X(t')=A c (t')[ Pll +A(t')}+eAl(t'), 



(31) 



and 



h+A(t)] 2 _[k + A(t)] 



■2eA c (t)[k- Pl \-p±] 



-P±- 

(32) 
We now solve the saddle point equations in terms of 
the ionization t' and rescattered t times for giving paral- 
lel p« and perpendicular (transverse) p± momenta. We 
focus on the solutions of the shortest orbits for returning 
electron considering our defined monochromatic field in 
Sec. MI Al The remaining sets of the orbits are strongly 
suppressed due to the wave-packet spreading. In par- 
ticular, we consider the cases when the electron tunnels 
around 7r/2 and rescatters at 2ir and when the electron 
tunnels at 37r/2 and rescatters at 37r. We refer to them 
as the first shortest pair (Pair s \) and the second shortest 
pair (Pair S 2), respectively. 

Figure [3] demonstrates the real part of the solution 
of saddle point equations of the rescattered electron, as 
function of the electron momentum pu , for several trans- 
verse momenta, for the Pair s \ . The upper and lower pan- 
els correspond to the rescattering and ionization times, 
respectively. The longer (dashed line) and the shorter 
(solid line) orbits, along which the electron returns, can 
be identified, since the real parts of t' and t are associated 
to the classical trajectories of an electron in a laser field. 
The long and short orbits coalesce for the minimum and 
maximum momenta p\\ for which the rescattering process 
described by saddle-point Eq. (|32|) has a classical coun- 
terpart. Beyond these momenta we are in the classically 
not allowed region, in which the yield decays exponen- 
tially. 



respond to e = 0, while panels (b) and (e) show e = 0.003 
and panels (c) and (f) depict e = 0.005. The upper and lower 
panels show the rescattering Re[t] and the tunneling Re[t'] 
times of the electron, respectively. The dashed and solid lines 
correspond to the long and the short orbits, respectively. 



For the homogeneous case, the classically allowed re- 
gion is centered around —2^JU\ ) (Fig. [5£a)). This center 
is determined by the most probable momentum the elec- 
tron may have at the time of the rescattering. In here, 
for vanishing transverse momentum p±, the cutoffs arc 



near 



'Up and 2 A 



. The classically allowed region 
shrinks and it results more localized around the center 
by increasing p±. For larger p± this region completely 
disappears. 

On the other hand, for the inhomogeneous case, the 
center of the classically allowed region is not anymore at 
around —2^/Up and it shifts towards different momentum 
values. Figure [5] demonstrates how the center and the ex- 
tension of the classically allowed region strongly depend 
on the value of e. For the case e = 0.003, this center 
is around 1.5-^/C/p and the cutoffs are near —\fU~p and 
5WUp for p± = (Fig. E^b)). Furthermore, for the case 
e = 0.005, this center is around 5y/Up and the cutoffs are 
near WUp and 6 y/Up~ for vanishing transverse momentum 
p± (Fig. EJc)). As overall, the classical allowed region 
decreases by increasing the inhomogeneity factor of the 
field e. In addition, this region shrinks much faster and 
more localized around the center as e increases. 

Figure [S] demonstrates the imaginary parts of the ion- 
ization t' (upper panels) and the rescattering t (lower 
panels) times, for the Pair s \. These solutions confirm 
the classical allowed and not allowed regions interpreta- 
tion. Indeed, they demonstrate that the imaginary part 
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FIG. 6: (color online) Imaginary part of the tunneling t' and 
rescattering t times for the electron, as a function of its paral- 
lel momentum pv, for the same parameters as in Fig. [5] Panels 
(a) and (d) correspond to e = 0, while panels (b) and (e) show 
e = 0.003 and panels (c) and (f) depict e = 0.005. The upper 
and lower panels show the rescattering Im[i'] and the ion- 
ization Im[i] times of the electron, respectively. The dashed 
and solid lines correspond to the long and the short orbits, 
respectively. 



FIG. 7: (color online) Real part of the tunneling t' and rescat- 
tering t times for the electron, as a function of its parallel 
momentum p«, when the electron tunnels out at around 37r/2 
and rescatters at about 37r, for the same parameters as in 
Fig. [5] Panels (a) and (d) correspond to e = 0, while pan- 
els (b) and (e) show e = 0.003 and panels (c) and (f) depict 
e = 0.005. The upper and lower panels show the ionization 
Re[t] and the rescattering Re[t'] times of the electron, respec- 
tively. The dashed and solid lines correspond to the long and 
the short orbits, respectively. 



of the rescattering times vanishes between the momenta 
for which the real part of the rescattering times coalesces. 
Beyond this region the imaginary t increases abruptly. 
This is a very clear indication that for the former and 
the latter regions the rescattering is classically allowed 
and forbidden, respectively. In addition, the Im[i'] ex- 
hibits a minimum near the center of the classically al- 
lowed region, even if there is no classically allowed re- 
gion. It means that rescattering is most probable for this 
specific momentum. On the other hand, the imaginary 
part of the ionization time Im[i'] of the electron is always 
non-vanished, due to the fact that the tunneling process 
has no classically counterpart. Furthermore, Im[£'] gives 
some indication on how the width of the potential barrier 
that the electron must overcome in the order to reach the 
continuum, varies with regards to e. In general, this po- 
tential barrier becomes wider as Im[i'] increases. In here, 
the Im[i'] decreases slightly as a function of e. 

Figure [Jj demonstrates the real part of the solution 
of saddle point equations of the rescattered electron, 
as function of the electron momentum p» , for several 
transverse momenta and for the Pair S 2- The upper and 
lower panels correspond to the rescattering and ioniza- 
tion times. 

For the homogeneous case, the classically allowed re- 
gion is now centered at a positive parallel momentum 
P\\ (Fig- EI a ) an d (b)) but it remains the same as the 
Pair s i (Fig. [5(a) and (b)). On the contrary, for the in- 
homogeneous case, the classical allowed region is not the 
same as the Pair s \. Indeed, the classical allowed region 
is now much larger than the case of the the Pair s \. In 



contrast to the previous case, the classical allowed region 
increases and more slowly shrinks and localizes around 
the center by increasing the inhomogeneity factor of the 
field e. For instance, for the case e = 0.003, the center 
of the classically allowed region is around ^\fUp and the 
cutoffs are near and 8.5 y/Up for p± — (Fig. [7(b)). 
In addition, for the case e = 0.005, this center is around 



G, 



'Up and the cutoffs are near 1.5 a 



'U„ and U\/U„ for 



vanishing transverse momentum p^ (Fig. [5(c)). 

Figure [5] demonstrates the imaginary part of times t' 
(upper panels) and t (lower panels) of the Pair S 2- These 
solutions confirm the classical allowed and not allowed 
regions interpretation of Fig. [7j In contrast to the previ- 
ous case, the Im[i] decreases slightly as a function of e, 
indicating the widening of the potential barrier through 
which the electron must tunnel out. 

In Figure [SJ we plot the real and imaginary parts of the 
solution of the saddle point equation of the direct electron 
(Eq. pop ), as a function of the electron momentum pii 
and for several transverse momenta. In all cases, the 
imaginary parts of tunneling time t' exhibit a minimum 
at the peak- field times cot — ir/2. This is due to the 
fact that the effective potential barrier through which 
the electron tunnels out is narrowest for these times. On 
the other hand, as the transverse momentum p±^ becomes 
larger the Im[t'] increases. This is consistent with the fact 
that the potential barrier widens in this case. 

For the homogeneous case, the imaginary time Im [£'] 
of the neighboring orbits are identical and has a mini- 
mum at pii = 0. In here, two neighboring orbits behave 
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the Im [£'] is different for each orbit. Indeed, this time 
is larger when electron tunnels at t = 2ir (dashed line) 
comparing to the case when electron tunnels at t = n/2 
(solid line). In addition, for the inhomogeneous case, 
the electron most probably reaches the continuum with 
a non zero momentum. Instead, the inhomogeneity fac- 
tor of the field e determines the most probable momen- 
tum. For instance, for e = 0.003 and e = 0.005, the 
electron will tunnels most probably with pii = 1.5 yUp 



and 



P\\ 



respectively. It means by increasing e 



the electron most probably will tunnel out with larger 
momentum. 



B. Spectra 



FIG. 8: (color online) Imaginary part of the tunneling t' and 
rescattering t times for the electron, as a function of its paral- 
lel momentum pv, for the same parameters as in Fig. [7] Panels 
(a) and (d) correspond to e = 0, while panels (b) and (e) show 
e = 0.003 and panels (c) and (f) depict e = 0.005. The upper 
and lower panels show the rescattering Im[i] and the ioniza- 
tion Im[t'] times of the electron, respectively. The dashed 
and solid lines correspond to the long and the short orbits, 
respectively. 
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FIG. 9: (color online) Tunneling time t' of the direct electron 
as a function of its parallel momentum py, for several trans- 
verse momenta p±. Panels (a) and (d) correspond to e = 0, 
while panels (b) and (e) show e = 0.003 and panels (c) and 
(f) depict e = 0.005. The upper and lower panels give the 
real Re[t'] and the imaginary Im[i'] parts of such times, re- 
spectively. The dashed and solid lines correspond to the long 
and the short orbits, respectively. 



symmetrically with respect to the laser electric field and 
a electron sees the same effective potential barrier for 
both orbits. Furthermore, the electron most probably 
will tunnel out with zero momentum. On the other hand, 
for the inhomogeneous case, the two orbits do not behave 
symmetrically with respect to the laser electric field and 



In this section, we compute the spectra of the direct 
and the rescattered ATI with Eqs. ([3]) and ((4]) using the 
saddle point method developed in Sec. Ill CI 

Figure [TU] presents ATI spectra for hydrogen atoms 
(I p = —0.5 a.u.) interacting with a monochromatic 
field of frequency uj = 0.057 a.u. and intensity / = 
3 x 10 14 W/cm 2 and for the cases with e = 0.0 (homo- 
geneous case), e = 0.003, and e = 0.005 as shown in pan- 
els (a), (b) and (c), respectively. In all panels the green 
line shows the direct ATI spectra. For the direct ATI, the 
cutoff is around the conventional 2U P value for the homo- 
geneous case and it extends to higher energy as a function 
of the inhomogeneity factor of the field e. The cutoffs of 
the rescattered ATI are in good agreement with the tra- 
jectories represented the previous Section (|III A 1|) . For 
the inhomogeneous case, the pairs in Fig. [Ha) lead to the 
largest cutoff (red lines) and the pairs in Fig. [T{b) give 
lowest cutoff (cyan). However, the latter case will give 
the final cutoff when we consider the total contributions, 
i.e. from all the pairs given in Fig.Q] as shown in dashed 
black line. In the spectra with total contributions, the 
pairs from Fig. QJb) contribute to develop more interfer- 
ences at the low energy part of the spectra. In conclusion, 
for the homogeneous case, the cutoff of the rescattered 
ATI is at the well known value of 10U P (Fig. rTOT a)). while 
these cutoffs are at around 17U P and 25U P for e = 0.003 
and e = 0.005, respectively. This fact is indeed consistent 
with quantum mechanical simulations |lO| . 



C. Angular momentum distributions 

In this section, we calculate the electron momentum 
distributions as a function of the electron momentum pn 
and p± to the laser field polarization and including both 
the direct and rescattered ionization processes. We as- 
sume that that the prefactors Eqs. (O, ^ and © are 
constant, in order to remove any momentum bias that 
may arise from such prefactors and having a clearer pic- 
ture of how the inhomogeneous strength influences such 
distributions. 



ionization rate increases. 
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FIG. 10: (color online) ATI spectra for hydrogen atoms [I p = 
0.5 a.u.) interacting with a monochromatic field of frequency 
u) — 0.057 a.u. (corresponding to a laser wavelength A = 
800 nm) and intensity 7 = 3 x 10 14 W/cm 2 for the case with 
e = 0.0 (panel (a)), e = 0.003 (panel (b)) and e = 0.005 
(panel (c)). The green colored spectra show the direct ATI 
and the rest show the rescattered ATI. The red and cyan 
colored spectra show the contributions from pairs given in 
Fig.rjja) and Fig.fljb), respectively. The dashed black colored 
spectra shows the total contributions from all pairs given in 

Fig.ni 



Figure [TT] shows these calculations for the direct ATI. 
For the homogeneous case (panel (a)), the distributions 
are centered at around vanished p\\ , while for the cases 



I 
with e = 0.003 and e = 0.005 are at around p 



and p 



This fact is indeed consistent with 



minima presented in lower panels of Fig. [9l For all cases, 
the overall momentum space populated by the electron 
is the same and the interference patterns are less pro- 
nounced as e increases. Another interesting feature is 
that the overall the yield increases with increasing of e, 
as shown in Fig. [TTj) . This is due to the fact that the 
Im[i] decreases as a function of e, and as a result, the 
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FIG. 11: (color online) Electron momentum distributions for 
direct ATI and for hydrogen atoms (I p = 0.5 a.u.) interacting 
with a monochromatic field of frequency u) = 0.057 a.u. and 
intensity I = 3 x 10 14 W/cm 2 . Panels (a), (b) and (c) show 
the case with e = 0.0, e — 0.003 and e = 0.005, respectively 

Finally, Fig. [12J shows the electron momentum distri- 
butions for the rescattered ATI. The top row shows the 
homogeneous case, i.e. e = 0, while the middle and bot- 
tom rows give the inhomogeneous cases with e = 0.003 
and e = 0.005, respectively. The left and middle columns 
represent the cases for Pair s \ and Pair S 2, respectively. 
The right column depicts when both pairs are consid- 
ered and added up coherently. For the Pair s \, the elec- 
tron momentum distributions are centered at around 




Up, p\\ = 2.0^/Up and p\\ = A.Q^/U V for 
0.003, and e = 0.005, respectively. For 
these centers are at around pii = 2.0y / [/p, 
P|| = 4.0^/Up and p\\ = Q.Q^/TT P for e = 0, e = 0.003, 
and e = 0.005, respectively. For homogeneous case, pan- 
els (a) and (d), both pairs have the same yield. On 
the other hand, for the inhomogeneous cases, Pair S 2 has 
much smaller yield in comparison to Pair s %. We also 
calculate the electron momentum distributions by con- 
sidering the contributions from both Pair s \ and Pair S 2- 
These calculations are shown in panels (g), (h) and (i) 
for e = 0, e = 0.003, and e = 0.005, respectively. In here, 
for the inhomogeneous cases, Pair s \ will dominate the 
total electron momentum contributions as it has higher 
yield. 



IV. CONCLUSIONS 

In this work, within the strong field approximation 
(SFA), we derive the transition amplitude for both di- 
rect and rescattered above-threshold-ionization (ATI) in 
an spatially inhomogeneous field. In particular, we show 
how the quantum orbits of the ATI (for both direct and 
rescattered processes) manifest themselves in spatially in- 
homogeneous fields. We show that in nonhomogeneous 
fields, the ionization and rescattering times are not the 
same for the neighboring cycles, even for a monochro- 
matic field. In here, the electron tunnels with different 
canonical momenta, which one leads to the lower energy 
cutoff and the other one develops a higher energy cutoff. 
We demonstrate that the energy cutoffs of the photo- 
electron spectra of the ATI extends to higher energies 
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as a function of the inhomogeneity strength of the field. 
For both direct and rescattered processes, the center of 
the electron momentum distributions shift toward differ- 
ent momentum depending on the inhomogeneous factor 
and the cycle considered. Furthermore, the overall yield 
of the direct ionization slightly increases as function of 
inhomogeneous factor. For the rescattered case, this is 
true if the electron tunnels at around 7r/2 and rescattered 
at about 2ir. For the case when the electron tunnels at 
3tt/2 and rescatters at 3n, however, the yield decreases 
as function of the inhomogeneous strength. 
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FIG. 12: (color online) Electron momentum distributions for 
rescattered ATI and for hydrogen atoms (I p = 0.5 a.u.) inter- 
acting with a monochromatic field of frequency ui — 0.057 a.u. 
and intensity / = 3 x 10 W/cm . The top, middle and bot- 
tom rows show the cases with e = 0, e = 0.003 and e = 0.005, 
respectively. The left and middle columns represent the case 
when the electron tunnels out around 7r/2 and rescatters at 
about 2-7T and when the electron tunnels at 37r/2 and rescat- 
ters at 37r, respectively, while the right column shows when 
we have both cases. 
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